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We offer a physicists' proof that center-vortex theory requires the area in the Wilson-loop area law 
to involve an extremal area. This means that unlike in string-theory-inspired models of confining 
flux tubes, area-law dynamics is determined by integrating over Wilson loops only, not over surface 
fluctuations for a fixed loop. Fluctuations leading to to perimeter-law corrections come from loop 
fluctuations as well as integration over finite-thickness center- vortex collective coordinates. In d = 3 
(or d = 2 -f 1) we exploit a contour form of the extremal area in isothermal which is similar to 
d = 2 (or d = 1 + 1) QCD in many respects, except that there are both quartic and quadratic 
terms in the action. One major result is that at large angular momentum ^ in d = 3 -I- 1 the 
center-vortex extremal-area picture yields a linear Regge trajectory with Regge slope-string tension 
product a'{0)KF of l/(27r), which is the canonical Veneziano/string value. In a curious effect 
traceable to retardation, the quark kinetic terms in the action vanish relative to area-law terms in 

1/2 

the large-^ limit, in which light-quark masses ~ Kp are negligible. This corresponds to string- 
theoretic expectations, even though we emphasize that the extremal-area law is not a string theory 
quantum-mechanically. We show how some quantum trajectory fluctuations as well as non-leading 
classical terms for finite mass yield corrections scaling with We compare to old semiclassical 

calculations of relativistic (massless) qq bound states at large £, which also yield asymptotically- 
linear Regge trajectories, finding agreement with a naive string picture (classically, not quantum- 
mechanically) and disagreement with an effective-propagator model. We show that contour forms 
of the area law can be expressed in terms of Abelian gauge potentials, and relate this to old work 
of Comtet. 

PACS numbers: ll.15.Tk, ll.15.Kc UCLA/03/TEP/11 



I. INTRODUCTION 



The most familiar and elementary concept of confinement in QCD is the area law for Wilson loops. For non- 
relativistic problems, such as the low-lying bound states of heavy quarks, the simple implementation of this idea is 
through a linearly-rising potential. But after thirty years of QCD it is still not at all clear how to address, in the 
continuum and for d ^ 2, the dynamics of light-mass qq bound states and similar relativistic problems where the 
notion of potential breaks down. In the early days of QCD, many authors (for example, ^) suggested that a good 
approximation in d = 4 would be to use an effective gluon propagator behaving like 

A,Ak') = {-9,. + ---)^^^ (1) 

for small momentum fc, where Kp is the string tension (the omitted terms in Eq. Q are longitudinal). The non- 
relativistic potential derived from this effective propagator is indeed linearly-rising^nd in the relativistic regime yields 
a linearly-rising Regge trajectory at large angular momentum £. (As discussed in the qq problem is relativistic for 
quark mass M if ^ Kpi.) But this effective low-momentum propagator cannot be entirely correct; for example, 
it cannot be used as an effective propagator between gluons themselves, which are screened and not confined. 

When phenomenological approaches such as p|, using Eq. were first taken up, there was little understanding of 
how confinement actually worked in QCD. Since then, many people have come to believe in the center-vortex picture 
of confinement. (For early work in the center- vortex picture, see [EQIEIEQl; an up-to-date review is in [^.) In this 
picture, center vortices (of co-dimension 2) can link topologically with a Wilson loop; with a spaghetti-like condensate 
of center vortices, the vacuum expectation value (VEV) over gauge fields of a Wilson loop results in an area law (for 
group representations where the center is non-trivially represented). The vortices have finite thickness, but as far as 
the area law is concerned, thickness effects represent perimeter corrections ignored in this paper. 

The purpose of the present paper is twofold: First, to argue, in physics language, that center-vortex dynamics leads 
to the prescription that, for a given Wilson loop, the area in the area law is extremal for that loop. Our second goal 
is to formulate the relativistic large-^ dynamics of quarks in this extremal-area center-vortex picture, and to compare 
this dynamics to earlier work, such as the effective-propagator model and various string models, in particular that 
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of We will find that the extremal- area picture agrees (to the extent that they can be compared) with but 
gives a different value for the Regge slope than does We will study mostly semi-classical dynamics, which should 
be appropriate for large angular momentum £, and indicate how integration over fluctuations of Wilson loops gives 
rise to corrections which vanish at large £. We will not address the full complications of relativistic qq dynamics, 
including quark spin, perturbative gluon exchange, and the like. Our focus will be on finding the asymptotic slope of 
the linearly-rising Regge trajectory which (we will show) arises in the extremal-area picture. 

Even before what we will call QCD-inspired models, such as that of 1], many authors had studied linearly-rising 
Regge trajectories in the context of the Veneziano model and various simple models of particles joined by strings; for 
example, see 0. The canonical relation between the string tension Kp and the Regge slope a'(0), based on making 
classical string theory agree with the Veneziano model, is: 

a'iO)Kp ^ ^. (2) 
zvr 

Other models give somewhat different values for the numerical constant on the right-hand side. More recently, string 
theorists have constructed models of d = 4 strings which could represent hypothetical strings joining quarks. 
There is no reason that the general QCD-inspired model should give the Veneziano/string theory Regge slope, since 
the Veneziano model, which pre-dates QCD, does not necessarily capture the true dynamical features of QCD, and 
there is a big difference between the properties of literal, physical strings which join quarks and the area-law dynamics 
proposed in the present paper, based on center vortices. Moreover, to the extent that string theory is conformally- 
invariant there is no room in the string-theory action for separate kinetic terms for quarks. (In an effective string theory, 
massive quarks can be represented as the ending of open strings on finitely-separated branes, yielding appropriate 
kinetic terms.) The difference between literal strings and center-vortex dynamics is that the area in the area law for 
literal strings can fluctuate even though the contour (Wilson loop) is fixed, and these fluctuations need to be taken 
into account, while for center- vortex dynamics the area in question is minimal, and there are no surface fluctuations 
for a fixed bounding contour. 

In Sec. ^ we discuss at some length the question of how one argues that the Euclidean center- vortex area law 
involves extremal areas. The standard naive argument leads to an area law, but it is far from clear what surface 
and area are to be used for a given Wilson loop. We adopt a center-vortex model which retains only those features 
necessary to discuss this question of the area law. The model consists of a condensate of infinitesimally-thin closed 
vortex lines {d = 3) or sheets (d = 4) a finite fraction of which have infinite length. These vortices are characterized 
by two lengths: A persistence length and a density of vortex pierce points per unit area. In this model we show 
that the area in the area law is indeed extremal. As stated above, finite-thickness vortices lead only to perimeter-law 
corrections. 

Once the extremal-area law is established, we face the problem (the Plateau problem) jlO| of finding the extremal 
area associated with a given bounding Wilson loop. In Sec. IIIII we remind the reader of some well-known obstacles 
to solving the Plateau problem. These include the behavior of extremal areas for poorly-behaved contours, such as 
space-filling contours, and the general existence of more than one extremal area for a given contour. We speculate 
that some of these peculiarities might actually be observable in QCD dynamics. For example, two parallel and equal 
circles separated axially have (at least) two extremal areas: The areas of the circles, and a catenoid joining them. One 
or the other is a global extremum, depending on the axial separation of the circles. The transition from one global 
minimum to the other may be related to the behavior of the adjoint potential, which is a breakable string. 

In Sec. IIVI we consider how to express a Euclidean extremal area in contour form. This is actually elementary 
but not commonly used (string theorists, for example, usually use standard integrals over the world sheet). The 
easiest case is d — 2, where any simple area is a global minimum, and can be expressed as a quadratic form in 
the contour. In higher dimensions the area law is simplest in Dirichlet form, which requires imposition of gauge 
conditions in the form of choosing isothermal coordinates. The Dirichlet form of the area is still quadratic in the 
variables describing the surface, but this is misleading, because this is so only in isothermal coordinates. We show 
that imposing isothermality leads to quartic terms in the surface variables, in any dimension greater than two. (In 
d = 2 isothermality is automatic.) After expressing the Dirichlet form of the extremal area as a contour functional, we 
take up the issue of how to enforce isothermality, and end up with a modified contour functional which is quartic in 
the contour variables. In certain cases these quartic terms may be dealt with as perturbations around a flat surface. 
This area contour functional, plus kinetic terms, describes a one-dimensional field theory of contours. This is to be 
contrasted to string-theoretic forms of confinement, which lead to two-dimensional field theories. It is far from easy 
(but not impossible; see ^2,]) to find the right form of such a theory for confinement, since standard bosonic string 
theory has a critical dimension of 26. 

In d = 3 isothermality can be identically satisfied by a choice of functions describing the surface, as Weierstrass 
and Enneper showed long ago. However, the Enneper- Weierstrass formulation does not make it easy to express an 
extremal area as a functional of its bounding contour, in general. As would be expected from the paragraph above, 
the Enneper- Weierstrass area form is quartic in its variables. 
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We also show in Sec. IIVI that the Dirichlet area law in contour form has an interpretation, in d dimensions, as 
a contour integral over d Abelian gauge potentials. Alternatively, inci = 3ord = 2 + l one may choose a single 
Abelian gauge potential and a complex scalar field as the degrees of freedom. These, as Comtet as shown, can be 
identified with Witten's T^l degrees of freedom for spherically-symmetric = 4 SU (2) gauge theory. We will pursue 
this connection further elsewhere. 

In Sec. 0we apply this formalism to some simple problems of (spinless) qq relativistic dynamics in rf = 2, = 1 + 1, 
and in Sec. I VII to the d = 3 + 1 problem of 0], (which is, at the classical level, actually a problem in c? = 2 + 1, since 
one spatial coordinate is inert). Given the action (including kinetic terms) as a functional of Wilson loops, one must 
integrate over all Wilson loops to construct Green's functions. An example of such an integral, for scalar "quarks" in 
Euclidean space, is the free quark propagator: 



where M is the scalar-quark mass; s,r are proper-time variables; and the integral (dz) is an integral over contours, 
subject to the constraint that these contours pass through a; at r = and y at t = s. For qq dynamics one introduces 
a quark and and anti-quark propagator such as in Eq. (Q, plus the appropriate area-law action, and integrates over 
both quark and anti-quark contours. As one expects, the d = 2ord=l + l dynamics is exactly equivalent to the 
well-studied and ancient problem of large- two-dimensional QCD, and our discussion of two dimensions is largely 
review. 

In this Section we take up the problem of circular qq orbits at large angular momentum i, using the extremal-area 
formalism, and compare to the results of ,lj, which used the effective propagator of Eq. and of 0, which used 
a naive string model. Very fortunately, it is quite elementary to find the extremal surface expressed in isothermal 
coordinates (a helicoid) and its area analytically, corresponding to the semiclassical orbits (a double helix) arising 
from solving the classical equations of motion. One surprising feature of the extremal-area model is that in the small- 
mass regime (meaning <C Kp£) the quark kinetic energy M7, where 7 is the usual relativistic factor, vanishes 
in comparison to the energy coming from the extremal-area action. This is rather different from the results of P], 
where M7 does not vanish relative to area-law energy. The significance of this vanishing is that, in some sense, the 
extremal-area law is rather like string theory which has no separate kinetic terms for quarks. 

Note that this small-mass regime does not require literal masslessness, which is in any case impossible in QCD, 

since even quarks of zero current mass undergo chiral symmetry breakdown driven by confinement itself, which leads 
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to mass generation with M ~ Kp . We will not attempt to model such effects here, and will in fact ignore area-law 
terms responsible for constituent-mass generation. But we do show that both non-leading classical terms and certain 
quantum corrections both yield additive corrections to the bound-state energy E which are of the same form as mass 

terms, that is, E = + 0{Ky^). 

As a result of the vanishing of quark kinetic energy relative to area energy for the extremal-area model, the Regge 
slope-string tension product of this model is precisely that of the Veneziano/string model of Eq. |2Jl. This slope 
differs significantly from that found in 0]; we find that this is because for any area- like model where quark kinetic 
energy is significant in the large-£ limit, the Regge slope is necessarily smaller than the Veneziano/string slope. We 
also show that identically the same classical area action governs the present work and that of These authors 
quote, but do not derive, the Veneziano/string slope as applicable to their model. The identity of classical actions 
does not mean that the string model of Q and the extremal-area model are physically identical, because presumably 
in the string model the world surface of the string is subject to the usual sort of quantum fluctuations that strings 
have, even when the Wilson loop it joins is fixed. But in the center-vortex picture of the present paper the nature of 
the fluctuations is quite different: For a given Wilson loop there is a unique surface, at least if the loop is sufficiently 
well-behaved, and perimeter-law fluctuations arise from the thickness of center vortices, not from surface fluctuations. 

In Section we also show how certain integrations over quantum fluctuations of the contours lead to corrections 
which vanish as £ — > cx). 



For simplicity, we restrict our explicit considerations to SU{2) gauge theory; this avoids complications with nexuses, 
baryonic Wilson loops, and similar topological elaborations occurring for SU{N), N > 3. Our center-vortex model 
consists of a vacuum condensate of infinitesimally-thick oriented vortices, of co-dimension 2 (so points in d = 2, closed 
lines in d = 3, closed two-surfaces in d = 4), characterized by two (presumably related) quantities: An areal density 
p of vortex pierce points per unit area and a persistence length I. The density p is the number of vortices penetrating 
an element of area, divided by that area, each pierce point being counted as +1. Each vortex can be considered as a 
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II. ARGUMENT FOR EXTREMAL AREA 
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closed random walk which avoids itself and other vortices (subject to the constraint that the density is p), and the 
persistence length is, as usual, the length over which the direction of the random walk decorrelates. 

As is well-known, the VEV of a Wilson loop is an average over vacuum gauge-field configurations of phase factors 
constructed from link numbers of the vortex configuration with the Wilson loop. For SU{2) the phase factor is 
exp(z7rL) for link number L, so that only the link number mod 2 counts. The total link number of a given vortex with 
a Wilson loop is found by adding numbers ±1 (depending on the relative orientation of vortex and Wilson loop) for 
each oriented intersection of the vortex with a surface spanning the Wilson loop. Because the linkage of an idealized 
zero-thickness center vortex and a Wilson loop is topological, it is tempting to say that the area-law part of the 
Wilson-loop VEV is stationary under changes of the surface spanning the loop when the loop is held fixed, and so 
the surface is extremal. This is, however, a seriously incomplete argument. Let us review the standard but sloppy 
argument for the area law in center vortex theory. We write the Wilson-loop VEV as: 



where Li = ±1 is the link number associated with the i pierce point of the spanning surface, and N is the total 
number of pierce points of vortices with the spanning surface. Some vortices intersect the spanning surface only once 
(or an odd number of times) and contribute -1 to the sum of terms in the loop VEV; others intersect not at all (or 
an even number of times) and contribute a phase factor -1-1, as if they were not present at all. 

Suppose that the Wilson loop under consideration is large, in the sense that all its spatial scales are large compared 
to the persistence length I and the density length p~^/^, and mathematically benign, for example, C°°, not self- 
intersecting, and so forth. Then the total number of pierce points N = pA of the spanning surface S, with area A, 
is large. Each pierce point is associated with a partial link number Li = ±1. A certain fraction ^ of these Li are 
independent and random; the remainder are link numbers of vortices which are small compared to the Wilson loop 
scales, and are not linked. One might then argue that the Wilson-loop VEV is found by assuming Gaussian statistics 
for the independent link numbers. (Actually, Poisson statistics, but this fine point is irrelevant to us). Then this VEV 
would be given by: 



with {Lf) — 1. Here Nind — S,pA is the total number of independent vortices. 

At this point an area law is apparently established, but it is far from clear what area is to be used. To say that 
the area A{S) in the area law is the area of an arbitrarily-chosen surface S makes no sense for the general spanning 
surface, for example, one whose scale lengths and area are astronomical in scale, although there is nothing wrong with 
using such a surface. The VEV is quite independent of which surface is chosen, contradicting the naive interpretation 
of Eq. ((Sj. Yet there must be some surface Swl, of area Awl, such that the VEV is given by this equation: 



This must be the area of some surface S spanning the Wilson loop; if S did not span it, some vortex linkages would 
not be counted. 

What went wrong is the assumption that all the Lj are actually independent, even for vortices very long compared 
to the Wilson loop size. Consider the vortex associated with a particular pierce point i. Let us consider a topological 
configuration consisting of two surfaces each spanning the Wilson loop, one of them being an arbitrary surface S that 
is far from the correct surface and the other the correct (but unknown) surface Swl ■ The surface S + Swl is a closed 
surface with no boundaries; we choose it to be oriented. Since all vortices are closed, if a vortex pierces the surface S 
it must pierce the surface S + Swl an odd number of times more. For simplicity, we state the argument explicitly 
for vortices which penetrate exactly once more. Consider the fate of a given vortex which enters the volume enclosed 
by -|- Swl by penetrating S. Upon exiting this volume it can penetrate either S or Swl again. If it penetrates S 
but not Swl it is not linked to the Wilson loop; otherwise it is. Evidently, the probability P that the given vortex 
will in fact penetrate the Wilson loop is P = A{Swl)IA{S). In effect, the mean-square link number (L^) is not 1, 
but P, and the area law does come out in the form of Eq. ©, and is independent of the originally-chosen surface S. 
Unfortunately, this argument does not tell us how to choose Swl- 

This argument holds unless the surface S is very close to the correct surface SwL^ in which case the link numbers 
and pierce points of the two surfaces are strongly correlated. Consider a new choice of spanning surface Swl + 5S 
which is very close to the correct surface, and calculate the Wilson-loop VEV with the surface Swl + 5S. The area of 
this spanning surface is Awl + 5 A. Either the link numbers Li of the correct surface Swl become those of the new 




(4) 



(e-^)=exp[— (i2)7V,„,] 



(5) 



{e^^^)=eM—{Ll)£.pAwL]. 



(6) 
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surface, without change, or new pierce points arise in pairs, with total hnk number 0. These new pierce points arise if 
in going from Swl to Swl + 5S an unhnked vortex which penetrated one of these surfaces twice does not penetrate 
the other one at all. But these effects cancel when averaged over many gauge configurations, and provided that the 
spatial scale of the deviation 5S is small compared to the persistence length, do not happen at all. The result is that 
the Wilson-loop VEV is precisely the same for the two surfaces: 

(e^-^) = eM~^ipA{SwL)] = exp[-yCp(^(^H/L) + 5A)] (7) 
and so 5A{Swl) = 0; the correct surface is extremal. 

III. EXTREMAL SURFACES AND CONFINEMENT: A HOPELESS PROBLEM? 

As mathematicians know, it is hopeless to give a completely general solution to the Plateau problem, because of 
the variety of contour and surface types that could be encountered. It is a theorem that any Jordan curve possesses 
at least one extremal area which spans the curve, but there can be infinitely many such areas, in pathological cases. 
There is no particular reason, for example, to suppose that the contour spanned by a minimal surface is at all regular. 
An example is a random walk, or something like Hilbert's space-filling curve. 

It could happen that such a curve generates an area which grows faster than R^, where i? is a maximal diameter 
of the contour. From Eq. 1221) below, one proves by elementary arguments the inequality for the minimal area A: 

A < ^RL (8) 

where L is the length of the contour and R a maximum diameter. For a contour generated by a random walk, 
L = R^a~^ where a characterizes the average step length for the walk. It is then possible, from Eq. (jHJ, that the 
area grows as fast as R^. However, if this rate of growth were actually achieved, it could not be important for the 
infrared behavior of Wilson loops, since the contributions of a loop averaged over gauge fields decays exponentially 
as exp{—KA) where K is the string tension. The mathematical reality is that we can only make progress with 
well-behaved contours, and we assume good behavior at every step of our arguments. 

Even with this assumption of well-behaved contours and surfaces, some interesting questions arise. There are many 
contours which have more than one extremal area. Consider, for example, two equal-size circles whose planes are 
parallel; depending on the separation between the circles (in units of the circles' radius) several extremal areas are 
possible, as shown in Fig. ^ We may think of the loops in Fig. ^as being idealized representations of nearby qq pairs 
as would occur in a hybrid meson. The mere geometrical existence of an extremal area does not mean it makes sense 
in a confining gauge theory. For example, the configuration of surfaces shown in Fig. is not sensible in SU{2) 
gauge theory, with the usual interpretation of the intersection of a center vortex (a closed line that, if linked K times 
to a fundamental Wilson loop, gives rise to the factor {—)^ for the Wilson loop). This intersection number K can be 
calculated by counting the number of signed intersections of the Wilson loop with an orientable surface spanning the 
loop. A simple center vortex can be linked once positively to both Wilson loops of Fig. I^K, and then it contributes 
a factor of unity to the product of Wilson loops, just as if it were linked to neither. But in Fig. |2K it is possible 
either for a vortex so linked to pass through the central plane surface and for it not to; in the former case, one would 
assign a factor of -1 to the product of Wilson loops, and in the latter, a factor of -1-1. This figure can make sense for 
SU{3) with the assignment of fluxes to surfaces consistent with the rules for associating a nexus world line with the 
line where three surfaces meet in the figure. 

Go back to SU{2); both Figs. C make sense. The same link numbers are found by counting penetrations of 
either spanning surface by a center vortex. But the areas to be used in the area law for Wilson loop VEVs are 
clearly different, and so apparently is the interpretation. For Fig. 03 the Wilson loops are correlated, in that the 
spanning surface depends on both loops, but in Fig. they are uncorrelated. Interestingly, as one would expect 
on physical grounds, when the separation of the loops is large enough (roughly equal to their common radius) the 
catenoid suggested by Fig. 03 is not a global minimum; only the independent planar surfaces of Fig. 03 are. This 
suggests the possibility of short-ranged van der Waals-like forces between two hadrons, not necessarily associated with 
meson exchange, and also has implications for the behavior of the screening law for adjoint Wilson loops, and area 
laws for iV-ality non-zero representations of SU (N) for N > A. We will explore these issues elsewhere. 

IV. FINDING THE CONTOUR FORM OF EXTREMAL AREAS 

Certain of our remarks apply in all dimensions, while some are specific to d = 3. Also of interest are elementary area 
formulas in d = 2, because their contour version resembles what we will find in higher dimensions to some extent, with 





A 
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FIG. 1: A. Two parallel circles with extremal surfaces meeting at the Steiner angle. B. Two circles joined by a catenoid. C. 
If the circles are far enough apart, the global minimum surface is the two planar surfaces of the circles. 



the contours acting as currents sources for a d = 2 Abelian gauge propagator. Later we will show that such Abelian 
gauge propagators arise in all dimensions. In another section we will review d = 2 dynamics, which is governed by a 
quadratic and hence soluble action. For d> 2 quartic terms arise from the isothermality conditions. 
The area of a surface is given by the usual expression: 



-I' 



where the <Ji, i = 1,2 are coordinates on the surface and r] is the determinant of the induced metric: 



Vij = 9al3 



(9) 



(10) 



Greek indices run from 1 to d. For now we will be concerned with a Euclidean embedding space, so g^p — Sa/s ■ We do 
not need to distinguish covariant from contravariant indices, but for notational clarity wc will always use contravariant 
notation for Greek indices and covariant notation for Latin indices. If wc impose the isothermal gauge conditions 



d(Ti dai 9(72 9(72 



A; 



dx" dx" 
dai 9(T2 



= 0. 



(11) 



the equations for an extremal surface of area A (defined as one with zero mean curvature everywhere) are derived 
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from a simple Dirichlet action: 



^D = l I cPa{d,x"f. (12) 



These equations of an extremal surface are V^a;"=0, which means that the coordinates can be extended to holomorphic 
functions of the complex variable a = Ui +102- We write this extension of real coordinates a;" to complex coordinates 

X" = .x" + iq°' (13) 

where x", are real and the X" are holomorphic functions of a. The Cauchy-Riemann relations relate the real and 
imaginary parts of the derivatives: 

d.x'^ = ey-^jg", a.g" = -eyS.x" (14) 

(with di = d/dcTi). By virtue of the Cauchy-Riemann relations the isothermality condition of Eq. (|11|) becomes: 

(d^X'^f = 0. (15) 

In Eq. H12|l we distinguish between the true area A and the Dirichlet area Ad , because in general, if the isothermal 
conditions are not imposed, the most one can say is that A < Ajj, as follows from: 

^ = \J dM[T.\9<rX'^\"^' -\T.(9.x-m'^'- (16) 



The Dirichlet area 



Ad^I I <fa\d,X^\\ (17) 



agrees with the true extremal area when isothermality is imposed. Note that one may write A^ 'm a. number of ways, 
for example. 



AD = \jd^Cj{d,q''f (18) 



as well as in the original form of Eq. 112|) . 

Two-dimensional area formulas illustrate these principles very easily. 

A. The area formula in two dimensions 

In d = 2 Euclidean space every area associated with a piecewise once-differentiable non-self-intersecting closed 
contour Xi{t) is minimal, so the required answer is just the usual area. 

We label the two coordinates Xi and choose ai = Xi (thereby fulfilling the isothermal conditions). In this case, Eq. 
I|22|l becomes the usual formula: 

A = ±i j dtcijXiXj (19) 

where the contour is given as a function of a parameter t, and the overdot indicates the derivative with respect to 
t. The area must be positive, of course, and the sign is chosen according to choices made about the direction of the 
contour. 

A second form, making clear the connection between c? = 2 QCD and an area law, is the d — 2 version of Eq. I|27|) 
and it states the area law as an action involving a massless vector propagator. This form is: 

A ^ (j) dxi (j) dx'jAij[x - x']-, Ay = (% H )(^) In |a; - a;'|. (20) 

(We allow for longitudinal terms, containing di or dj and indicated by the dots; they contribute nothing to the area.) 
Of course. Ay is the d — 2 massless vector propagator, so that large-iV d = 2 QCD is truly an area law, as is 
well-known. 

One can show directly that Eq. (|20|l gives the usual area by using Stokes' theorem on each contour integral, plus 
the fact that the Laplacian on Ay gives a delta function, to arrive at ^ = / d^a. 
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B. Contour formulas in higher dimensions 

The contour form of the Dirichlet area is found by Stokes' theorem and the observation that 

=e,,a,(g"a,a;"). (21) 

One contour form of yl^i is then: 

ADc = \i da\"d,x'' (22) 



2jr 

where F is the contour spanned by the minimal surface. We give it a new name Ad^ for reasons which will become 
clear later. 

The functions q" are not immediately known in terms of the a;", but they can be eliminated by using the Cauchy- 
Riemann relations in the standard integral for an antiderivative G: 

G{a) = — I da'G'ia') ln{cr' - a) + const. (23) 
2vri /r 

where G" = dG/da' . The integral gives zero if a lies outside the contour F; the unspecified constant depends on how 
the contour is constructed. 
We write 

ln((T' - cr) = In |cr' - cr| + 17: 7 = arctan[^|— ^] (24) 
and observe the Cauchy-Riemann relations: 

ei-jdj In \a' - cr| = d^. (25) 

It is then straightforward to integrate by parts twice in Eq. (|23|l . keeping track of endpoint terms which can survive 
if the contour crosses the logarithmic branch line, and find: 

G = Gi+iG2^^ <j> d(T[ In \a' - a\e,jdjGi + ^ <j> da[ In \a' - a\diGi + const. (26) 

This equation allows us to express the , which are the unphysical imaginary parts of the surface coordinates and 
correspond to G2 in Eq. H2ti|l . as integrals over derivatives of the physical coordinates (corresponding to Gi). By 
expressing in terms of dix" through equations of the type of Eq. H26() one finds that Auc is given by: 

Aa^^^j dcrAx'^ia) <j> da'^djx'^ia') In |cr' - a\. (27) 

Eq. (|27|l is apparently a contour-integral form for the Dirichlet area. However, more work is needed, because in Eq. 
(|27|l one can specify independently functions a;" (a) and a contour on the world sheet Oi = (Ti{t), where the parameter 
t labels points along the contour. The contour in target space is a;"(t) = x'^[ai{t)]. For a fixed target-space contour 
the area Adc can take on a continuum of values. Insuring that the contour area Adc agrees with the Dirichlet area 
A]j requires that isothermality be imposed. 

C. Dealing with isothermality 

We now show that if the isothermal conditions are imposed that the contour form of the area, Aoc, depends on the 
contour ai{t) only to the extent that the original form of the Dirichlet area (that is, Ajj of Eq. (|17|l ) depends on the 
contour, when the functional dependence of x" is held fixed. This means that for a variation 6ai of the a contour, 
the target-space contour changes by: 

5x"it) = 5a,{t)d,x"{a,). (28) 
The contour variation of the Dirichlet area A]j is: 



SAd ^ <p e^jda.SajidkX^y. (29) 
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As for Adc one finds, by using the Cauchy-Riemann relations, its variation: 

6Anc = '2 j) dai6akeudix°'dkx°'. (30) 

The isothermal conditions are: 

dix'^dkx" = ^6ik{dpx"f (31) 

and one finds that dAoc = <5^r>, given isothermality. Note that the form of 6Auc^ when isothermality is imposed, 
yields reparameterization invariance, since by definition reparameterization means that 

eijdfJiSdj = 0. (32) 

Of course, isothermality is imposed only along a contour. But since the isothermal relations are analytic in character, 
they extend from the contour to its interior. In order that the contour integral for Adc be useful, we must impose 
isothermality on the functions appearing in it. A useful starting point in d = 3 is the Enneper-Weierstrass construction, 
which reduces isothermality to an identity. 

In d = 3 the isothermal gauge conditions are imposed, as Enneper and Weierstrass showed long ago, by writing 

d„X = {l-f)g- d„Y = -i{l + f)g- d„Z^2fg (33) 

where g and f^g are meromorphic in the domain of a containing the minimal-surface contour (so / can have poles if 
they are cancelled by zeroes of g). 

The function g changes under conformal transformations, and so for a wide class of contours, those conformally 
related to a standard contour, one can choose g — 1. Alternatively, for this class one can choose a standard contour 
in the a plane, for example, the unit circle, and determine an appropriate g. 

Note that two-dimensional areas are a trivial sub-case of the Weierstrass construction. For simplicity, choose the 
x^-x'^ plane as the d = 2 space; then in the Weierstrass construction, take f — 0, g — 1. The anti-derivatives giving 
2^1, 2^2 yield (up to a translation of the origin) Xi = ai, and the area formula of Eq. H34() is in the non-contour form 
/ d^x. This is, as already discussed, readily translated into a contour form. Although in d = 2 one chooses / = it 
is not necessary to take g — 1. But ii g ^ 1 the specified contour in the xi^X2 plane is not the same as the contour 
in the cti, cr2 plane. 

There are, of course, infinitely many contours which can be associated with an extremal area described by given 
functions f,g, provided that the contours lie within the domain of harmonicity of these functions. For any of these 
choices, the extremal area is given by: 

A = j d^a\g\\l + m\ (34) 

Weierstrass long ago showed 'll^ how, in principle, to solve the extremal-area problem for a simple non-intersecting 
closed contour composed of straight-line segments. However, the solution has a character rather like the closely-related 
problem of Schwarz-Christoffel transformations, in which it is not at all easy to relate the parameters of the solution 
to the parameters of the contour. 

Deviations from two-dimensionality might be profitably described with /, g written as Taylor series. Even a simple 
case such as keeping 17=1 and saving only a linear term in / (quadratic and cubic terms in the coordinates) results in 
a rather complicated surface known as Enneper's surface (see Fig. O. 

An alternative to using the Enneper-Weierstrass functions /, g, is expansion of the coordinates X" — (AT, y, Z) 
themselves as Taylor series and then imposing isothermality. We choose the origin of coordinates so that the surface 
passes through the origin at cr = and write: 

^ = E>";^ = E^-";^ = E>"- (35) 
1 1 1 

One easily calculates the isothermality condition of Eq. in terms of these coefficients: 

N 

= '^{ajaK + PjPk + 1j1k)\j+k=n+i (36) 
1 
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FIG. 2: A portion of the Enneper surface. 



It is straightforward to see that the iV = 1 term in Z merely corresponds to a tilt of the surface, and can be set to zero. 
Correspondingly, we set a2, /?2 — 0. We can always choose a\ — if3i — 1, as we did in d = 2. Under the assumption 
that all the 7jv are of order of a smallness parameter e, there is a perturbative construction of isothcrmality: 

"3 - i/?3 = -^72; "4 - iPi = -7273; "5 - iPs = -7274 - ^("3 + (^3+ Ts)' ■ ■ • (37) 

One more condition can be imposed; the simplest is to set the Enneper- Weierstrass function g to unity, so that 
aN + t-Pn = for iV > 2. The holomorphic function / is just: 

/ - \daZ. (38) 

Now we have just enough degrees of freedom for three contours, namely, two tr-space contours Uiit) and the contour 
defined by ReZ(o'i(t)), whhich can be chosen freely. It is far from trivial, in general, to find explicitly the target-space 
contours {t) , although this can be done perturbatively. One may say that this difficulty is the real difficulty of the 
Plateau problem, and we have nothing new to say about it. But from the point of view of QCD dynamics, what is 
required is not to associate a set of target-space contours to areas, but to integrate over these contours, or areas, in a 
way consistent with external boundary conditions. Integrating over the a contours and the z contour will serve just 
as well, provided that no difficulties are encountered in imposing the boundary conditions. 
Let us define new functions i^i, F2 as the real and imaginary parts of an antiderivative of 



(39) 
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and then calculate the area from Eq. H27() . We specify a unique antiderivative by choosing F{a = 0) = 0. We assume 
that / begins with a linear term and so F begins with a quadratic term in a (a linear term in F simply amounts to 
a tilt). 

The area Adc, which we need no longer distinguish from the true extremal area A, is: 

A = jda,[U,,a, + '-{Z*d,Z - Zd,Z*) + '-{F*d,F - Fd,F*)]. (40) 

As advertised, there are d = 2 constructions of a type already encountered in Eq. H19|l or equivalently Eq. H2()|l . 
plus a quartic term. Note that the action (area plus kinetic terms) connected with these degrees of freedom is not 
a d = 2 theory, but a c? = 1 field theory of the contours, although these contours are coupled by a massless d — 2 
propagator, as in Eq. H20() . 

The contour form of the Dirichlet area can be interpreted in terms of Abelian gauge potentials, as we now take up. 



D. Area and Abelian gauge potentials 

In any dimension d > 2 the area can be expressed as a sum of d contour integrals over d Abelian gauge potentials, 
which are derived from the complex coordinates X". Moreover, we can implement gauge transformations on these 
potentials and show, as expected, that they do not affect the area. 

In view of the area relations of Eq. H12II18|) . we can express the area as an arbitrary linear combination of an x-form 
and a (/-form and integrate by parts: 

A^\<j> da,€ij [Xx°'djX°' + (1 - X)q"djq°'] (41) 

where A is a numerical parameter. (A different A can be used for each dimension, but we need not indicate that 
explicitly.) By using the Cauchy-Riemann relations one changes this to 

A = (42) 

where 

Af = ir^.x" - x^d^q"] + (i - ^)a,(x"g") (no sum) (43) 

One sees that a choice of A is equivalent to choosing a gauge for the d gauge potentials A"; we will work now in the 
gauge A = 1/2. In this gauge we find that: 

A^ = liiX'^ydiX" - X^idiX"")*] (no sum). (44) 
8 

Another way of writing these A = 1/2-gauge potentials is: 

Af ^le.jdilX'^l^ {no sum). (45) 



Straightforward calculation shows that 



Al +A^ = ie,,a, + [F*d,F - Fd,F*] (46) 



and, according to Eq. (|44|l . 



Al='-[Z*d,Z-Zd,Z*]. (47) 

The area constructed from the sum of the three gauge potentials is just that of Eq. 1)40(1 . 

It is interesting to note that we can reinterpret the degrees of freedom in the area law of Eq. H40|) as consisting 
of one Abelian gauge potential coupled to a complex scalar field Z, with quadratic and quartic terms in d^Z {Z 
can be expressed in terms of daZ by Eq. 1(231) ). These are just the degrees of freedom in Witten's description 
of spherically-symmetric S'[/)2) gauge potentials, which he used to describe certain multi-instanton configurations. 
Comtet [l3| has shown that the equations for extremal surfaces in = 2 + 1 can be reduced to Witten's instanton 
equations. We will discuss this connection in detail elsewhere. 
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V. TWO-DIMENSIONAL DYNAMICS 

We begin with a description of the area law in Euchdean {d — 2) space, and then, to set the stage for higher 
dimensions, complete the problem by considering d = 1 + 1. 

A. d = 2 

For two spinless quenched quarks, described by a field (f>{x), the Euclidean two-point mcsonic Green's function is 
written: 

G(m-w) = (O|0^0(w)0V(")|O) = / ds ds' {dx){dy)e-^ (48) 

Jo Jo J 

where {dx){dy) is a path integral over quark contours Xi(r) and anti-quark contours Ujir') and the action I is: 



2 Jo ^'^^'^Y io '^^'yj^Y^^^^"'^~ Jo '^^'^^j^^^J + — dT'€,jy,yj. (49) 

In writing this action we have omitted self-mass terms, coming from quadratic contour integrals involving two quark 
(or two anti-quark) orbits. These are necessary for insuring gauge invariance, and so must implicitly be understood 
to be included when questions of gauge invariance are brought up. In Eq. Kr is the string tension, M is the 
quark mass, and the boundary conditions on the path integral are: 

Xi{T = 0) = Mi; Xi{T = s) = Vi\ yj{T' = 0) = vj; yj{T' = s') = uj. (50) 

The mesonic spectrum is revealed by the poles in momentum of the Fourier transform of this Green's function. 
The two propagators in the Green's function are represented by the proper-time integrals. In two-dimensional QCD 
Kn = g^Cfi, where is the coupling and Cr is the quadratic Casimir eigenvalue. 

Evidently the action is quadratic in the field variables, and thus exactly solvable. However, one must be careful 
about the sign of the area term, since contours of both orientations occur. This Euclidean action is simply the 
Minkowski action for a particle in a constant magnetic field, and the Green's function and particle spectrum are found 
from the classical solutions to the equations of motion. The M = case can only be studied by postponing the limit 
until all calculations have been performed. 

The classical solution is conveniently written in terms of the complex variables 

X = xi + ix2] U — ui+ iu2\ V = vi+ iv2 (51) 

and is: 

-iKrtr/M _ i 

^-^ = (^-^) ,-.KWM_i - (52) 



The action for, say, the quark (proper time s) has the kinetic part 



, Kf^s iv^-u,y Ms 
AM'I-C 2 



where in this and future equations we will use the shorthand 

The anti-quark contributes an identical kinetic term, with s — > s'. The area part of the action for the quark is: 

— dre^jX^Xj ^ Kr[-{viU2-V2Ui) + - S)\. (55) 

This represents the sum of two areas. The V1U2 — U1V2 part is the area of a triangle whose vertices are the origin, Ui, 
and Vi , and the remainder is the area bounded by the quark contour going from Ui to Vi and a straight line from Vi to 
Ui. The anti-quark area term has the same form with two changes: 1) the sign of the V1U2 — V2U1 (triangle) term is 
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FIG. 3: Geometric picture of the areas associated with the quark and the antiquark; each is bounded by a portion of a circle 
and a straight hne between the initial and final positions of the qq pair. 

The contributions of the straight lines cancel in the action, so that the action depends on the contour of the overall area. 



reversed; 2) s is replaced by s' in the remainder. The sum of the quark and anti-quark terms is just the area bounded 
by the sum of the quark and anti-quark contours (see Fig. EJ. The total action, then, is of the form: 

'^(4M^^3^ + — + 4(1-C) + (56) 

As required, this is translationally- invariant. 

The area terms (last terms on the right in Eq. (|5(j|l ') have an interesting geometric interpretation. Elementary 
geometry shows that the quark area is that of a circle of radius 

2\sil{KRs/2M)y ^^^^ 

less an area bounded by a portion of the boundary of the circle (the term linear in s) and a chord drawn from Ui to 
(the term -S). This is shown in Fig. 01 for the case tt < Krs/M < 2tt. For Krs/M in the range {2t:N, 2'k{N + 1) 
the circle is covered N times. It is, of course, this infinitely multiple covering of areas which allows the buildup of 
poles in the Fourier transform of the Green's function. 

In forming the mesonic Green's function by integrating over proper times, one must be careful, because the sign 
of the area term changes according to whether the quark is to the left or to the right of the anti-quark (as viewed 
with respect to the straight line joining u and v). This necessitates the introduction of absolute- value signs at various 
places, and complicates the analysis. The same issues occur in d = 1 + 1 dimensions, which we take up now. 



B. d = 1 + 1 



First we must give the translation from Euclidean to Minkowskian variables. In so doing, we must distinguish 
covariant from contravariant indices in the usual way, and we will use Greek indices for Minkowski-space variables. 
The change from Euclidean to Minkowskian space in higher dimensions is entirely analogous to the d = 1 + 1 discussion 
that we give next. 

Since in the Euclidean version it was possible to identify the two target-space coordinates with the two surface 
parameters, this must also be possible in Minkowski space; consequently, the surface coordinates are Minkowskian, 
described as a^,a^ with metric (cr")^ — (cr^)^ = ct^. Evidently the Euclidean area formula is to be replaced by the 
Minkowski-space formula 



A™ = / dV(-77)i/2 (58) 
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where ij is the determinant of the Minkowski-space world sheet metric, of signature (1,-1)- This "area" is real, as 
it must be. Of course, just as in d 2 we may identify the world sheet with the target space, and the area can be 
written in contour form as: 

A,n= (f da°'eapxxd^x^. (59) 



Manipulations similar to those for d — 2 show that the complete Minkowskian action /,„ for the qq system is, analogous 
to the Euclidean action of Eq. H49(l . 

- /m = Y J dTXc^x"' + Y J '^'''ypy'^ + + s') + Kr j dx°' j dyc,A{x - y). (60) 
The qq Green's function is constructed from the integral: 



oo 



m{x)m)) = / ds ds'e^'-. (61) 
Jo Jo 

In Eq. H6Ul the quantity A is an inverse of the d'Alembertian, and must be real in order to have a real area law in the 
action. We cannot use the Fcynman propagator Ap, which has both real and imaginary parts: 

AFix-y)^^\n\xf^leix'). (62) 

Instead we use only the real part, which is also the old Feynman- Wheeler propagator, half the sum of the retarded 
and advanced propagators. So we keep (as in m in d = 3 + 1) only the 0- function, and identify 

A{x - y) ^ ^e[{x - y)']. (63) 

Actually, as points out, this is gauge-equivalent to replacing the 6'-function by —0[—{x — y)^]. 

We need not give further details of this well-known c? = 1 -t- 1 problem which is, of course, the standard hyperbolic 
problem of acceleration in a constant electric field of strength Kr, so that sine and cosine in the d — 2 formulas 
becomes sinh and cosh. The energy E of the quark (or anti-quark) is constant, and given by: 

E = Mj + Kr\x' -y^\ (64) 

where the coordinate x^{y^) is the quark (anti-quark) spatial variable, and 7 is the usual relativistic factor for either 
the quark or the anti-quark (which always have equal and opposite velocities). The absolute value sign on x^ — y^ is 
the Minkowski analog of the signs we mentioned for the Euclidean case, and occurs because the field strength changes 
sign as the quark and anti-quark pass each other. WKB quantization leads to the standard formula i?^ = ttK^N for 
large integral N (in which limit the mass M can be neglected). 



VI. THREE-DIMENSIONAL DYNAMICS 



First we study the d = 3 + 1 problem of qq bound states at large angular momentum I, which at the classical level 
is actually a d = 2 -|- 1 problem because the orbits are planar, and compare it to two earlier approaches to the same 
problem. Then we give a scaling law for corrections to the classical result arising from certain quantum fluctuations. 

We are interested, as in the previous section, in the bound-state masses of qq systems considered classically. Of 
particular interest is the limit of zero quark mass M, where we will find a novel effect in which the quark kinetic 
terms in the action vanish in the limit. In actuality, even massless quarks protected by chiral symmetry acquire a 
constituent mass, since confinement breaks chiral symmetry, but we will ignore these effects here to make our main 
point, which concerns the asymptotic slope of the QCD linearly-rising Regge trajectory. We will, as in the d = 1 -I- 1 
case, ignore explicit mass-generation terms involving two quark or two anti-quark contours, but these are important 
for maintaining gauge invariance. 

Classical orbits are sensible for large angular momentum which is all that we will study here. In this limit, the qq 
orbits lie in a spatial plane (and, of course, evolve in time), with the result that one spatial coordinate can be ignored 
and so a classical problem in d = 3 -I- 1 can be formulated in d = 2 -I- 1 Minkowski space. We expect massless quarks 

1 /2 

to get constituent masses of order KJ. , and this effect leads to corrections to our results which vanish at large £. 
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There are several works from long ago that have also considered this problem, and we will be concerned with 
comparing the results found here to two of them P, . In both these two works and in the present investigation 
the bound-state masses lie on a linearly-rising Regge trajectory in the massless quark limit. It turns out that in this 
limit the Regge trajectory has slope q;'(0) obeying the usual Veneziano/string relation a'{0)KF ~ l/(2vr) for both the 
extremal-area picture and for Ref. 9], which is based on a naive string picture. However, the third model JJ, which 
is based on the effective propagator of Eq. has a different and smaller value for this product. 

At the classical level, at least, the present work is virtually identical in its predictions to that of Ref. In fact, the 
extremal-area action is identical to an action (not explicitly quoted in 9]) which yields the results of that reference. 
But the present work and any form of string theory must differ at the quantum level, as we have discussed earlier. 



A. The effective-propagator model 

Consider first the effective-propagator model of Ref. Ij. In this model the action corresponding to Eq. H60|l is: 

-/™ = Y^'rfri„i" + dT'y0f + ^{s + s') + ^l^drl^ dr'i"(r)y„(T'M-(x - yf] (65) 

where only the real part of the propagator has been saved. (We study here only the case of equal mass M for the 
quark and anti-quark; the general case can be found in 1].) 

The contours chosen in 1] to solve the classical equations of motion in c? = 3 -t- 1 for the effective-propagator 
dynamics, with action as in Eq. H65|l. are: 

= ( — , pcos( — ), psin(^), 0); (66) 
s s s 

V = ( — , -pcos(^),- psm( — ), 0). 
s s s 

Here p is the common radius of the quark or anti-quark orbit, and parameters are chosen so that at r = s one has 
= (t, pcost/), psint/), 0) (and similarly for the anti-quark). For the purposes of finding bound-state energies, 
boundary conditions at t = are irrelevant. The z-coordinate is fixed, and so the problem is really in d = 2 + 1. 
From now on we omit writing the inert z coordinate. 

As explained in to calculate the total energy E of and angular momentum £ of a bound state one can set 
s = s', t = t' , 4> = (f)' . These prescriptions are consistent with the requirement that we are working in the bound- 
state center of mass, so that the energy E is actually the mass of the bound state. We then impose four conditions 
corresponding to the four external variables of the problem: s, t, p, (f). Following the imposition of these conditions we 
introduce a derived quantity to, obeying <f> — ojt where uj is the constant angular velocity of the quark or anti-quark. 
The common quark linear velocity is v — pui. 

First, a steepest-descent integration over these proper time variables in the Green's function shows that s = 
t/7, s' = t'/7 where 7 = (1 — v^)~^^'^ is the usual relativistic factor. (We will see that only the kinetic terms depend 
on the proper times.) The second and third conditions follow from the definitions of energy and angular momentum 
as derivatives of the action: 

The fourth condition comes from the equations of motion. There is an azimuthal component of these equations which 
is identically satisfied because the propagator we use is half the sum of the advanced and retarded propagators. The 
radial component of the equations of motion is perhaps most easily found from the equation 

op 

0=^1. (68) 

which is a way of saying that p — 0. From these four conditions one eliminates s, M7, and ui to find the energy as a 
function of £; this is, not unexpectedly, a linearly-rising Regge trajectory. 

For the effective-propagator model the action is (ignoring terms which vanish as t 00): 

M - 0^(1? , M,t2_ 202 r 2. 2, . />,M 

" " " s '^'^^ 7^ + ^ + + p^sinmi (69) 
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Here f is the positive root of the equation 

i2^2 = 2p2(l + cos(C0)). (70) 

The meaning of ^ is that it is the value of |r — r'|/s at which the 0- function in the effective-propagator action of Eq. 
(|74|l vanishes. A httle algebra shows that one may write ^ = 9/ujt where 9 is the positive root of 

e"^ =2v^{l + cos9). (71) 

For massless quark and anti-quark, with u — 1, one finds 1\ 9 ^ 1.478, and in the non-relativistic limit of small w, 
9 = 2v. 

It remains to carry out the action derivatives in Eqs. I|67|) . In so doing, the actions is expressed as a function of 
s,p,(j),t. After the derivatives are taken, one may then substitute s = t/^,(j) — cut. We are only interested in the 
rclativistic (massless) case, so where it is safe to do so we set v — 1 in the resulting equations. This may not be done 
in the product M7 until this combination has been eliminated among the three equations. We find: 

E = 2Af7 + — Ci; £ = + ^^2; M-ftu = KpCs (72) 

where the Ci are the values of functions Ci{v) at 1; = 1; these values depend on 9{v = 1). A couple of lines of algebra 
yields a linearly-rising Regge trajectory: 

E'-C,K,l;C,^^^§±^. (73) 



B. The extremal-area model 



For the extremal-area action, the only difference is that the area-law interaction term Ia has a ^-function with a 
different argument as prescribed in Eq. (|27|) (when transcribed to Minkowski space): 

-Ia^^ f dr r dT'x'^{T)y^{T')9[-{a - a')'] (74) 
^ Jo Jo 

where cr(r) ((t(t')) is the quark (anti-quark) contour on the world sheet, and x"(r) = x"[cr(r)] is the target-space 
quark contour. In the non-relativistic limit (where 9 = 2v for the effective-propagator model) the extremal-area model 
and the effective-propagator model agree. In the massless limit where puj = v = I, and using = ujt, the area part of 
the action is: 

- ImA\v=l = (75) 

which is of the expected form. 

Note that the orbits for the effective-propagator model (Eq. H66|l ^ are helices, and taken together they form a double 
helix. These helices form the boundaries of a well-known extremal area, the helicoid, and it is not surprising that the 
extremal-area model has a solution involving this particular extremal surface. The helicoid has an axis running along 
the proper-time direction. The contour z"((T'^,tT^) spanned by the minimal surface can be very simply expressed in 
isothermal coordinates: 

z =( , -sm( )cos( ), -sm( ) sm( )). (76) 

S S S S S 

It is elementary to check that z"((t*', ct^) is harmonic and isothermal. 

The next step is to relate the Wilson-loop contour as expressed in world-sheet variables to the kinematic variables. 
For the quark orbit we identify the proper time r with a'"\ and for the anti-quark we identify r' with cr^. The other 
world-sheet variable (j^ is independent of proper time, and takes on values differing by a sign for the quark and 
anti-quark. We have: 



^0 



a" = T, ^ = sin(-!^ — ) (quark orbit); cr" = t', ^ = — sin( ) (anti — quark orbit). (77) 



The resulting helical orbits are essentially those of Eq. H66f) used for the effective-propagator action. 
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The real difference between the effective-propagator model and the extremal-area model now shows up. Because 
p4>lt = V, where v is the linear velocity of the quark or anti-quark, we can write for the boundary values of on the 
contour: 

(ha^ IT 
M^O: — = iarcsint; -> ±-. (78) 

One now sees, by consulting the argument of the ^-function in the extremal-area action of Eq. (|74|l and comparing 
it to the 0- function of the effective-propagator action of Eq. (|65|l . that the effective limiting value of the angle 9 
introduced in Eq. H71|) is n at zero quark mass, or unit quark velocity. This angle 9, as discussed in Q, measures 
the retardation (or advance) of signals exchanged between quark and anti-quark. The rate of approach to the limit is 
such that 

. „ const. , ^. 

sin 6* ^ (79) 

7 

(we do not need to know the constant). Some straightforward algebra shows that at 9 — n the equation of motion 
is different in structure from that of the effective- propagator model, as given in the third equation of Eq. (|72|) . since 
the coefficient of Kp vanishes proportional to sin^?. Given the rate of approach of vanishing, the equation of motion 
for the extremal-area model is: 

M-fu = (80) 
7 

for some (irrelevant) constant C3. This is easy to solve: 

M -> : Mj^ (MpRpC'^y^^ -> 0. (81) 

The kinetic terms drop out of the action, leaving only the area term, in the massless limit. 

Of course, if the massless limit were in fact attainable, and in this limit the quark kinetic energy vanished, classical 
calculations would be highly suspect. An infrared massless particle does not possess the kind of localization expected 

1 /2 

for classical physics to be valid. In fact, a constituent mass of 0{Kp ) will be generated, so there is some localization. 
We should require that M"/p > 1 as a minimum requirement of localizability. This means, using Eq. (|81|l and the 
easily-derived scaling p ^ (i/Kpy^^ , which holds in the massless limit, that: 

M.MMp7.,)"»>i;,= M^f">^ (82, 

1 /2 

which, with M ~ Kp , simply means large £. So we will continue to use the massless limit as a formal device, 
expecting corrections from real masses which vanish at large £. 
This extremal-area action is: 

For comparison with the work of Callan et al Q we note that this action has the integral representation: 

- I^A = 2Kf r dRie - {R4>fY'^. (84) 
Jo 



C. Numerical values for the three models 



The constants for the effective-propagator case were worked out long ago l] and yield 

E"^ 9.91A>^ (effective propagator). (85) 

This yields a Regge slope-string tension product a'{0)Kp of about 0.101, rather small compared to the 
Veneziano/string theory product of l/(27r)=0.159. 

For the extremal-area model in leading order in I/7, we find Ci — t:, C2 = 7r/2, and to all intents and purposes 
C3 = 0, because of the extra factor of I/7 in Eq. (|80|l . The values of Ci, C2 are the same as found in the elementary 
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string model of Callan et al'wl. (This reference does not quote a value for C3.) Mathematically, at least, this is no 
coincidence; differentiation of the action for the extremal-area model as given in Eq. H84|l yields precisely the integrals 
for E,£ of Ref. But the physical interpretation is presumably quite different, since a string theory differs at the 
quantum level from the extremal-area theory. The values of the Ci of the extremal-area model yield, from Eq. I|73|l . 
the usual Veneziano/string relation: 

E'^ = 2'KKFi (extremal area). (86) 

One might expect that corrections to these leading-order results are of 0(7""'^) ~ in view of Eq. 1)81(1 and 

the leading-order scaling p ^ [1/ KpY^'^ . In fact, because of the specific form of the area-law action in Eq. H83|l . 
determined in part by dimensional considerations, it turns out that the leading corrections are 0(7^^) ~ 
leading to: 

E ^ {2^KFtf'^[l + 0{j^^^)]; (87) 

that is, the leading correction is just a mass term. It appears that at least some quantum corrections coming from 
orbital fluctuation in the (classically irrelevant) z-direction have the same scaling, as we discuss next. 



D. Scaling laws at large angular momentum from quantum fluctuations 

We estimate contributions from z-fluctuations away from planar circular orbits, as expresssed by the quartic terms 
in the area of Eq. (|4U|I . These estimates will be done in Euclidean space, which should have no effect on simple scaling 
laws. The final result identifies a smallness parameter e related to these deviations from planarity: 

e ^ (iff p2)-i/2 ^ ^-i/2_ (88) 

We will take into account only the 72 term in the Taylor series expansion for Z (Eq. I(35|l ) to estimate the corrections. 
Keeping only the 72 term corresponds to F = — 7|cr'^/12 in Eq. (|39|l . From Eq. (|40|1 for the area, multiplied by Kp 
to find the area part of the action, one finds that the quartic contribution £4 to the action scales as: 

/4-XfP^|72|^ (89) 

The quadratic contributions to the action are all positive, and so the functional integral over 72 leads to an estimate 
of the typical value of I72I: 

I72I - {Kfp')-''\ (90) 

Then Z ~ pe^^^ and the area part of the action (which is dominant over kinetic terms) coming from the F terms in 
Eq. (gni scales as: 

h Kpz^ Kpp^t. (91) 

To relate this scaling to the angular momentum ^, we have already seen that the mass Mg of a relativistic bound state 
with large £ scales as: 

Mt ^ {Kpif^ ^ Kfp, (92) 

and that: 

p-{^y/';er^r'/\ (93) 
Finally, the correction to the bound-state mass behaves like: 

Mi - (Kpiy^^ii + oi£-^/^] = {Kpey^^ + o{Ky^). (94) 

More accurate estimates will be given elsewhere. 
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VII. SUMMARY 

In this paper we argued in some detail for expressing the dynamical area law of QCD as a law of extremal surfaces, 
based on the topological mechanism of confinement in center- vortex models. In principle the area of an extremal 
surface can be expressed as a functional of the contour (Wilson loop) spanned by the surface. Confined-particle 
dynamics and bound states follow from functional integrals over these contours, which amounts to a d = 1 field 
theory. This is unlike earlier attempts to deal with area-law dynamics, in which string surfaces associated with 
tubes of confined fiux have d = 2 fiuctuations which must be integrated over, even for a fixed Wilson loop. These 
fluctuations, existing over the whole world surface of the flux tube, lead to severe complications 0. But in the 
center- vortex picture, fiuctuations have effect only on subdominant (perimeter) terms and are much more easily dealt 
with. 

We also pointed out some of the challenges associated with this law, such as the existence of several extremal areas 
for certain contours, and discussed some questions of the behavior of extremal areas for the sort of non-differentiable 
contours likely to occur in the functional integral over Wilson-loop contours. Certain more complicated contour 
configurations may have to do with adjoint string breaking and related phenomena. 

The contour form expression for an extremal area has an interpretation in terms of Wilson loops of d Abelian gauge 
potentials, in d dimensions. There is a natural definition of a gauge transformation, leaving the area unchanged, on 
each of these gauge potentials. It is interesting to note that in d = 3 the degrees of freedom can be reinterpreted as one 
Abelian gauge potential plus a complex scalar field, which happen to be the degrees of freedom identified by Witten 
[T^ in his spherical ansatz for SU{2) gauge fields. Comtet ^3 has shown that the SU{2) self-duality equations for 
these degrees of freedom are precisely those extremalize a surface in d = 2 + 1. What the connection is between our 
work and this earlier work remains to be seen; it is under investigation. 

These contour formulas are simple generalizations of area formulas in d = 2, which express an area as the interaction 
of a massless d — 2 gauge propagator with currents running along the Wilson-loop contour, just as required by d = 2 
QCD. However, their simplicity is deceiving, since the contour form of the area law holds only for the Dirichlet integral, 
which expresses a true extremal area only under the imposition of isothermal gauge conditions. Isothermality leads 
to expressions both quadratic and quartic in the contours, and so the extremal-surface form is an interacting field 
theory of these d = 1 contours. 

We gave one simple example in d = 3 -f 1 (actually applicable to d = 3 -|- 1 problems) of a classical extremal-area 
law, in which the quark and anti-quark contours were helices and the extremal area a helicoid. We compared the 
large-£ results of this center- vortex model with two earlier models, one an effective- propagator model and one a 
naive string model and found that the present results differ from those of but agree with those of even 
to the extent that the actions of ^ and of the present paper are mathematically identical. However, this does not 
mean that the physical content of these two models are the same; presumably at the quantum level the string model 
would involve functional integrations over string fluctuations, while in the center-vortex picture such fluctuations are 
missing. 

In the heavy-quark limit, all approaches including that given here agree with one another and with the conventional 
quantum mechanics of a linearly-rising potential. But there is a curious effect in the limit of quarks whose mass M 
approaches zero. One generally expects that in this limit the quark kinetic energy M7, where 7 is the usual relativistic 
factor, is finite. This is, for example, what happens in the effective-propagator model. But in the extremal-area 
model, retardation effects lead to an extra factor of 1 in the force term of the qq equation of motion, which has 
the consequence that A-Ij ^ M^/^ in the limit, and thus vanishes. The quark kinetic terms in the action no longer 
contribute, and only the area-law part is left. This, of course, is what string theory would prescribe, and in fact the 
extremal-area model studied here yields precisely the Veneziano/string result a'{0)KF — l/(27r). Generally, in any 
model where quark kinetic terms survive the massless limit, this product is smaller than l/(27r). 

Under the circumstances where the semi-classical developments of this paper are applicable, such as for bound 
states of large angular momentum i, there is a smallness parameter ^^^/^ which allows one to expand deviations from 
three-dimensional contours, dominant at large £, into the fourth spatial dimension. We made some estimates which 
showed that such fluctuations yield corrections vanishing at large £, as expected. 

The dynamics of extremal surfaces has just begun to be explored, and future work will probe more deeply into 
several of these matters. 
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